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SUMMARY 

Maize is one of the most important agricultural crops in Croatia, and was selected for research of the 
effect of climate warming on yields. The Decision Support System for the Agrotechnology Transfer 
model (DSSAT) is one of the most utilized crop-weather models in the world, and was used in this 
paper for the investigation of maize growth and production in the present and future climate. The 
impact of present climate on maize yield was studied using DSSAT 4.0 with meteorological data from 
the Zagreb-Maksimir station covering the period 1949-2004. Pedological, physiological and genetic 
data from a 1999 field maize experiment at the same location were added. The location is 
representative of the continental climate in central Croatia. The linear trends of model outputs and the 
non-parametric Mann-Kendall test indicate that the beginning of silking has advanced significantly 
by 1-4 days/decade since the mid-1990s, and maturity by 4-5 days/decade. It also shows a decrease in 
biomass by 122 kg/ha and in maize yield by 216 kg/ha in 10 years. 

Estimates of the sensitivity of maize growth and yield in future climates were made by changing the 
initial weather and C0 2 conditions of the DSSAT 4.0 model according to the different climatic 
scenarios for Croatia at the end of the 21st century. Changed climate suggests increases in global 
solar radiation, minimal temperature and maximal temperature (x 1 -07, 2 and 4 °C, respectively), but 
a decrease in the amount of precipitation (x0-92), compared with weather data from the period 
1949-2004. The reduction of maize yield was caused by the increase in minimal and maximal 
temperature and the decrease in precipitation amount, related to the present climate, is 6, 12 and 3%, 
respectively. A doubling of C0 2 concentration stimulates leaf assimilation, but maize yield is only 1% 
higher, while global solar radiation growth by 7% increases evapotranspiration by 3%. Simultaneous 
application of all these climate changes suggested that the maize growth period would shorten by 
c. 1 month and maize yield would decrease by 9%, with the main reason for maize yield reduction in 
Croatia being due to extremely warm conditions in the future climate. 



INTRODUCTION 

Climate changes have different intensity in different 
regions and so their impact on agricultural production 
varies. Thus, there is a need to research the impact of 
climate change on agriculture at national and regional 
levels. The impact of weather on crop growth, 
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development and yield can be best represented by 
crop-weather models, which facilitate the study of the 
relationship between weather, climate and crop yield. 
These models can be empirical-statistical or dynamic- 
deterministic. The empirical-statistical models use 
regression methods, while the dynamic-deterministic 
ones require morphological and physiological plant 
observations and meteorological measurements. Cur- 
rently, it is the dynamic-deterministic models that are 
mostly used for estimating yield and productivity, 
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leading to better control and management of agro- 
ecological systems, and also to better understanding of 
the cause-effect relationship in the soil-plant-atmos- 
phere system. The main application of crop models 
is in the study of the impact of climate change on 
crop production, whereas the operational application 
focuses on crop yield forecasting (Eitzinger et al. 
2008,2009). 

The aim of the U.S.A. IBSNAT project 
(International Benchmark Sites Network for Agro- 
technology Transfer), which started in 1982, was to 
develop a computer simulation method that would 
be widely applicable in farming, environmental 
protection and economics (Tsuji & Balas 1993). The 
Decision Support System for Agrotechnology Trans- 
fer (DSSAT) project (Hoogenboom et al. 1995) was 
developed within IBSNAT, which resulted in the 
DSSAT computer programme for PCs being made 
available to as many users as possible. The DSSAT 
programme for different crop types connects crop 
type, soil and weather conditions in a single pro- 
cedure. It includes the main crop types intended for 
human consumption, divided into three groups: 
cereals and maize, leguminous plants and root and 
tuber crops. Each crop group has its own basic 
simulation model, which is then adapted to a 
particular crop. The most widely used are the 
simulation models for maize and wheat, under the 
common name of CERES (Crop-Environment Re- 
source Synthesis). The CERES models for maize 
and wheat are some of the most frequently applied 
crop-weather models in the world, not just in Europe, 
and this justified its choice in this work. 

Maize and winter wheat are the most important 
agricultural crops in Croatia. Agricultural land 
amounts to 0-56 of the Croatian landmass and in the 
period 1998-2007, maize was produced on 0-32 of the 
arable area, winter wheat on 0-21 and other cereals on 
0-13 (Simanovic et al. 2008). 

The maize vegetation period in Croatia spans the 
warm season, from May to October; therefore, unlike, 
for example, winter wheat, maize growth and yield are 
not directly affected by snow cover. Most studies on 
the impact of climate change on maize production 
have dealt with agronomy but this paper will 
emphasize research from the meteorological aspect. 
Crop-weather models are not widely used in Croatia 
and this paper will help raise awareness of the im- 
portance of such applications for a better understand- 
ing of the effect of climate changes on agriculture. 

MATERIALS AND METHODS 

Meteorological, maize phenological and yield data 

The meteorological data used by the DSSAT model 
are daily values of maximal and minimal air tempera- 
ture, the amount of precipitation and the global 




Fig. 1 . Location of Zagreb City in Croatia. 



solar radiation. The field experiment was carried out 
at the farm of the Zagreb University Faculty of 
Agriculture in 1999, and the meteorological data used 
in the analysis were taken from the nearest meteo- 
rological station, Zagreb-Maksimir (45°49'N, 16°2'E, 
123masl), located at about 650 m from the field 
experiment site (Fig. 1). For the sensitivity analysis of 
weather conditions on maize production in the present 
climate, long-term time series of daily meteorological 
data were used for the period 1949-2004. The global 
solar radiation measurement at Zagreb-Maksimir 
covered too short a period for the necessary analysis, 
and data were taken from the Zagreb-Gric station 
(45°49'N, 15°59'E, 157masl) for the period 1949- 
1986 and the Zagreb-Horvatovac station (45°50'N, 
16°0'E, 182masl) for the period 1987-2004. These 
meteorological data are representative of the con- 
tinental climate in central Croatia, with a mean 
annual air temperature of 10-6 °C and a precipitation 
amount of 856 mm. 

For the validation of the DSSAT model, long-term 
time series of maize phenological stages were used 
from the Bozjakovina phenological station (45°49'N, 
16°18'E, HOmasl) for the period 1980-2004. 
Bozjakovina is about 20 km to the east of the 
Zagreb-Maksimir station. 

Historical data for the entire Croatian area for the 
period 1885-2008 were used to estimate the variability 
of total maize yield and maize-planted area in Croatia 
(Stipetic 1991; Simanovic et al. 2008). 

Soil data 

In accordance with IBSNAT recommendations 
(1990a, b), samples of the vertical pedological profile 
were taken for the chemical and physical analysis 
of the soil about 2 weeks before sowing, which was on 
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3 May 1999. The analyses were performed at the 
Pedology Department of the Faculty of Agriculture in 
Zagreb. The physical analysis by pedological profile 
layers determined the required pedological data: water 
content, field capacity, maximal saturation content, 
moist bulk density and wilting point. Chemical 
analyses gave the soil pH values per layer, the pro- 
portion of humus, ammonia, organic and total 
carbon, the ratio of nitrate in moist and dry soil and 
the mass of phosphorus and potassium. According to 
the current Croatian soil classification, the soil used to 
carry out the experiment belongs to the enteric brown 
soil type. Medium silt loam was determined by 
mechanical soil analysis. All measured soil data are 
presented in Vucetic (2006). 

Field experiment 

The field experiment at the farm of the Faculty 
of Agriculture in Zagreb lasted from 3 May to 
16 October 1999. One of the aims was to study how 
different weather conditions from year to year 
changed the maize grain yields despite constant 
practice management and soil conditions. The 1999 
field experiment was managed following typical local 
practices, with only one 'treatment' and without tile 
drainage or irrigation. In Croatia, maize is not 
irrigated and soil is watered only by precipitation. 
Therefore, it was assumed that in the DSSAT model 
only meteorological conditions differed from year to 
year in the period 1949-2004, whereas all other initial 
conditions remained constant. Thus, temporal maize 
yield variability over different weather conditions has 
been assessed. 

In the autumn and spring 1999 field experiments, 
248 kg N/ha was applied to the soil. The autumn 
fertilization (25 November 1998) was applied via deep 
furrow (0-25-0-30 m) and consisted of 450 kg of NPK 
fertilizer and 350 kg of urea. Spring fertilization 
was carried out on 15 April 1999 (300 kg of NPK) 
and 25 May 1999 (100 kg of limestone ammonium 
nitrate (KAN)). The maize cultivar PIO 3901 was 
planted to enable the results to be compared with 
research results in neighbouring countries. The size 
of the area sown was 8-5 x 20 m, which allowed for 
12 rows of maize at 0-7 m row spacing. Plants were 
0-19 m apart within each row and plant density was 
80000 plants/ha. Sowing depth was 50 mm. As per 
IBSNAT requirements (1990a) three biomass harvests 
and a final harvest were carried out. The biomass 
consists of the stem, leaves, leaf sheath, ear and husk, 
while yield is the grain yield. The determination of the 
moist matter mass, the drying of particular plant 
parts, the measuring of their dry matter and the 
chemical analyses were performed at the Special Plant 
Production Department of the Faculty of Agriculture 
in Zagreb. All measured crop data are presented in 
Vucetic (2006). 



Input data for the DSSAT model also include 
genetic coefficients based on the maize phenological 
phases. PI is the growing degree-day (8 °C tempera- 
ture threshold) from seedling emergence to the end of 
the juvenile phase, P2 is the photoperiod sensitivity 
coefficient, P5 is the growing degree-day (8 °C 
temperature threshold) from silking to physiological 
maturity, G2 is the maximum kernel number per plant 
and G3 is the potential kernel growth rate (mg/kernel/ 
day). 

DSSAT model 

The DSSAT model for maize (Jones & Kiniry 1986; 
Ritchie et al. 1990; Hoogenboom 2000; Jones et ah 
2003), besides simulating maize growth, development 
and yield also assesses the commencement of the 
phenological phases that depend on genetic coeffi- 
cients and weather conditions, leaf, stem and root 
growth, the development and growth of vegetative 
and reproductive organs depending on the phenologi- 
cal phases, soil water balance and water use by crop, 
soil nitrogen transformation, uptake by the crop and 
partitioning among plant parts. By including nitrogen 
and water balance in the model, it is possible to 
optimally use fertilizers and irrigation (where avail- 
able) to realize nutrition and water storage in the 
plant. Thus, in the last version of the DSSAT 4.0 
Cropping System Model (CSM) the different modules 
individually simulate the various processes that affect 
the land unit including plant growth, soil processes, 
soil-plant-atmosphere interface and management 
practice (Jones et al. 2003). 

In this paper, the DSSAT 4.0 model has been used 
for the estimation of maize production in the present 
and changed climate. In order to investigate the 
sensitivity of weather conditions on maize growth 
and yield in the present climate, the pedological 
conditions, maize genetic conditions and soil manage- 
ment were assumed to be the same as in the 1999 
typical field experiment. Thus, in the DSSAT 4.0 
model calculation, the Zagreb meteorological data 
were changed year by year during the 1949-2004 
period. This 56-year time series was used to estimate 
the beginning of silking and maturity dates, grain 
yield, kernel mass, aboveground biomass, maximum 
leaf area index (LAI), grain N and total N uptake. As 
the DSSAT 4.0 model includes a module of water 
balance, output data also included evapotranspira- 
tion, evaporation, soil evaporation and runoff. Then, 
a linear trend analysis was applied to the output 
results to determine the climate and maize production 
variability in Croatia in the present climate, in the 
period 1949-2004. 

Statistical analysis 

One of the methods for estimating the existence of a 
significant trend is the non-parametric Mann-Kendall 
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Fig. 2. (a) Long-term time series of total maize production (10 6 t) and area planted to maize (10 6 ha) and (b) mean maize 
yield (t/ha) as well as linear trends and 5-year running averages of total and mean maize yield in Croatia during the period 
1885-2008. Linear trends are significant at P<005. 



rank test, based on the values of individual elements 
of the series and the position of these elements in the 
series (Mitchell et al. 1966). If a linear trend exists, the 
values should increase or decrease chronologically. 
For significant linear trends (P<0-05), Sneyers (1990) 
progressive analysis was used for determining the 
beginning of a linear trend and its significance. 

A linear regression analysis of predicted maize 
phenological stages was done on the observed 
phenological data. For both statistical analyses, the 
STATISTICA 6.0 statistical package was used. 

Climate scenarios 

In order to investigate the sensitivity of maize 
production in future climate, the initial weather and 
C0 2 conditions of the DSSAT model were modified 
using the results of the different climatic scenarios for 
Croatia at the end of the 21st century (Vucetic 2008). 
The climate change scenarios were prepared by the 
pattern-scaling technique using the following global 
climate models (GCMs): ECHAM4/OPYC3 from the 
Max-Planck Institute for Meteorology, HadCM3 
from the UK Hadley Centre Climate Prediction and 
Research and CSIRO-Mk2 from the Australian 
Commonwealth Scientific and Industrial Research 
Organisation. A validation analysis in the Czech Re- 
public (Dubrovsky et al. 2005) showed these GCM- 
based standardized scenarios to be a good choice for 
representative climate changes scenarios in the con- 
tinental climate. Therefore, these GCMs were used for 
central Croatia. The mean global temperature (T a ) 
values were estimated by the one-dimensional climate 
model MAGICC for a middle-emission scenario 
(SRES-B2 + Al), combined with intermediate climate 
sensitivity (change in r G = 2-5°C; IPCC-TGCIA 
1999,2007). Thus, the middle increments are based 
on an average of increments related to a combination 



of middle climate sensitivity and SRES-B2 emission 
scenario, and middle climate sensitivity and SRES-A1 
emission scenario. 

RESULTS 

Trend analysis of long-term maize land and yield data 

Maize production has always been very important in 
Croatian agriculture. The total arable maize area 
increased from 350000 ha in 1885 to 550000 ha in the 
mid- 1930s, but the total production was very low 
(Fig. 2a). Yield was relatively constant up to about 
1945 but started to rise in the early 1950s and this 
growth continued until the beginning of the last war 
in Croatia (1991-95), when both maize growing area 
and yield suddenly dropped by 40-50%. Currently, 
the total arable maize area amounts to c. 300000 ha 
with a 2 Mt grain produced. The linear trend in 
10-year periods indicates a significant increase in grain 
yield (by 0-8 t/ha) for the period 1949-2008, which is 
faster than in the period 1885-2008 (Fig. 2b). 

Analysis of meteorological and maize phenological data 

The mean monthly temperature regime from March 
to October 1999 was characterized by a warmer period 
than the 1961-90 norm; in particular, it was very 
warm in March and extremely warm in September 
(Table 1; Katusin 2000). The positive deviation of 
mean monthly temperature in 1999 was mainly 
associated with a mean minimal temperature above 
normal. The previously mentioned months were also 
dry. May was rainy, but the summer precipitation 
amount was average. 

During maize growth and development, daily 
weather variations are more important than monthly 
variations; the daily cycles of precipitation and 
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Table 1. Mean monthly and annual global solar radiation ( SR), maximal ( t max ) and minimal temperature ( t m i„) 
and precipitation amount (Pr) for Zagreb-Maksimir during the period 1949-2004 and in 1999. Also given is a 
comparison of monthly and annual mean air temperature (t) and precipitation amount (Pr) in central Croatia in 
1999, compared with 1961-90, according to Katusin (2000). Temperature: C, cold; N, normal; W, warm; VW, 
very warm; EW, extremely warm. Precipitation: D, dry; N, normal; R, rainy 
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Fig. 3. Daily cycles of precipitation amount (Pr, mm), maximal (< milx , °C) and minimal temperature (/^n, °C) for 
Zagreb-Maksimir during the field maize experiment at the farm of the Faculty of Agriculture at Zagreb from May to October 
1999. 1, sowing (3 May); 2, silking (15 Jul); 3, maturity (15 Sep); 4, harvest (16 Oct). 



maximal and minimal temperatures during the 1999 
field experiment are displayed in Fig. 3. 

Rainfall amounted to 253 mm from sowing (0 on 
the BBCH code of Lancashire et al. 1991; Meier 
1997, 3 May) to silking (65 on BBCH code, 15 July), 
23 mm above the 30-year average and to 150 mm 
from silking to physiological maturity (87 on BBCH 
code when black dot/layer is visible on the base of 
kernels, 15 September) 9 mm below the 30-year 
average. In the sensitive reproductive growth stage, 
during the second part of July 1999, there was 
enough rain (61 -6 mm) and the mean minimal and 
maximal temperatures were 16-2 and 25-9 °C, respect- 
ively. The maximal daily temperature was above 



30 °C only on 20 and 21 July 1999. There were 
mainly optimal weather conditions for maize fertiliza- 
tion and the beginning of the grain filling phase. 
The warmest period was between 6 and 10 August 
1999 when the daily maximal temperatures were 
greater than 30 °C and the maximum temperature of 
33-7 °C occurred on 10 August 1999. After that, 
54 mm rain alleviated these unfavourable temperature 
conditions. 

A comparison of the phenological stages from the 
1999 field experiment with the 24-year normal for 
Bozjakovina shows an earlier beginning of milk ripe 
and wax ripe in 1999 (Table 2). Silking occurred on 
the same date as the long-term average. 
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Table 2. Linear regression coefficients of monthly and annual mean (t), maximal ( t max ) and minimal temperature 
(t, ni „), and precipitation amount (Pr) and maximal precipitation (Pr maxx ) for Zagreb-Maksimir during the 
period 1949-2004. Significant linear trends (P^O-05) are in bold 



Month 


t (°C/decade) 


'max (°C/decade) 


'min (°C/decade) 


Pr (mm/decade) 
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Table 3. Mean dates of maize phenological stages ( MEAN) and standard deviation ( s.d.) for Boijakovina in the 
period 1980-2004 and for 1999, and for the field experiment ( FE) at the farm of the Zagreb Faculty of Agriculture 
in 1999. The BBCH-code ( Biologische Bundesanstalt, Bundessortenamt and Chemical Industry ) is a decimal code 
system used to identify the phenological development stages of a plant ( Lancashire et al. 1991; Meier 1997) 



BBCH 


Phenological stages 


MEAN (date) 


s.d. (day) 


1999 


FE 


0 


Sowing 


29 Apr 


4 


28 Apr 


3 May 


10 


Emergence 


12 May 


4 


11 May 


15 May 


65 


First silk emerged 


15 Jul 


7 


9 Jul 


15 Jul 


75 


Milk ripe 


13 Aug 


10 


6 Aug 


5 Aug 


85 


Waxy ripe 


31 Aug 


12 


19 Aug 


20 Aug 


87 


Physiological maturity 


19 Sep 


16 


10 Sep 


15 Sep 


120 


Harvest 


3 Oct 


12 


16 Sep 


16 Oct 



It was important that the 1999 field experiment 
experienced 'typical' weather justified its use for 
investigating the impact of present and future climate 
changes. In order to establish the climate variations 
in present climate, the linear trends of monthly 
temperature and precipitation were defined for the 
Zagreb-Maksimir station for the period 1949-2004 
(Table 3). The linear trend analysis and Mann- 
Kendall text indicate that the greatest significant 
increase in mean monthly temperature in August 
(0-95 °C/decade) is mainly the result of a significant 
increase in minimal temperature. In general, a signi- 
ficant positive trend of monthly minimal temperature 
occurred from March to October and in January and 
of maximal temperature only in March and May. The 
monthly minimal temperature increased faster than 
the maximal and there was a significant annual 
increase in minimal temperature (0-4 °C/decade). 
There was no significant linear trend in monthly or 
annual precipitation in the period 1949-2004. 



Model calibration and validation 

The calibration consisted of determining the five 
genetic coefficients to adjust phenological stage timing 
and yield parameters. After each run of the DSSAT 
4.0 model, the coefficients controlling silking and 
maturity dates were readjusted automatically. The 
following genetic coefficients have been determined 
for the PIO3901 hybrid: P1=200°C; P2 = 0-10 and 
P5 = 740°C; G2 = 625 kernel/ear and G3 = llmg/ 
kernel/day. 

The variables used for validation were silking 
and maturity dates, grain yield, kernel mass, maxi- 
mum LAI, biomass, harvest index, N grain and total 
N uptake. The predicted values of the DSSAT 4.0 
model, derived from the meteorological, pedological, 
physiological and cultivation data measured during 
the 1999 field experiment are presented in Table 4. 
The performance validation of the DSSAT 4.0 was 
also done through a linear correlation between a 
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Table 4. Predicted values according to the DSSAT 4.0 
models and the observed values collected during the 
field experiment and from the pedological profile at the 
farm of the Faculty of Agriculture at Zagreb in 1999 



Description 


Observed 


Predicted 


Predicted/ 
Observed 


Sowing date: 3 May 1999 






Silking (day of 


196 


197 


101 


year) 








Maturity (day of 


258 


254 


0-98 


year) 








Grain yield (t/ha) 


131 


11-9 


0-90 


Kernel mass (g) 


0-347 


0-483 


1-39 


Maximum LAI 


4-8 


4-6 


0-96 


(m 2 /m 2 ) 








Biomass (t/ha) 


22-4 


21-8 


0-97 


N grain (kg/ha) 


138 


180 


1-30 


N uptake (kg/ha) 


177 


202 


114 


Harvest index 


0-58 


0-54 


0-93 



long-term time series of predicted maize silking and 
maturity dates, and the observed phenological data of 
maize silking and maturity dates in the Zagreb area in 
the period 1980-2004 (Fig. 4). 

Predicted mean values of maize components 

According to the theoretical calculations achieved 
by the DSSAT 4.0 model for the period 1949-2004, 
silking starts, on average, on 19 July, with a standard 
deviation of 6 days (Table 5). The beginning of silking 
and maturity were exceptionally early in 1950, 2000 
and 2003, which were extremely warm years in 
Croatia. 

Physiological maturity starts usually 70 days after 
silking with an 18-day deviation from average. The 
mean grain yield is about 12-4 t/ha and the average 
biomass is 22-2 t/ha. Both these predicted values can 
deviate by c. 1 -5 t/ha. In 1952, 2001 and 2003, the yield 
was under 9-0 t/ha. 

In the Zagreb area, the average kernel mass has 
been estimated at 0-46 g and the average LAI at 
4-4 m 2 /m 2 . The DSSAT 4.0 model simulated the mean 
amount of nitrogen in harvest seeds to be 183 kg/ha 
and the total nitrogen uptake during the vegetation 
period to be 257 kg/ha. 

Linear trends of maize yield components in present 
climate 

To determine the relationships of major changes in 
maize growth, development and yield over the 56-year 
period with weather conditions, the linear trends of 
particular predicted maize components were analysed. 
Silking started significantly earlier (— 1 -4 days/ 



decade), as did the achievement of maturity (—4-5 
days/decade; Fig. 5 and Table 6). There was a 
significant linear decrease in maize yield (—216 kg/ 
ha/decade). There was a slight decrease in kernel mass 
(— 0-01 g/decade), aboveground biomass (—122 kg/ 
decade) and harvest index (— 0-01/decade). There was 
no significant trend in kernel number/ear and maxi- 
mum LAI. The increase in total nitrogen uptake in 
maize by 0-6 kg/ha and the decrease in nitrogen in 
harvested seeds by — 0-03 kg/ha/decade is not signifi- 
cant in the Zagreb area (Table 6). A positive trend 
in evapotranspiration and soil evaporation (c. 3 mm/ 
decade) and in runoff (0-6 mm/decade) has been 
noticed during the vegetation period. An increase in 
runoff could cause a significant increase in soil 
erosion, which is mainly associated with the cultiva- 
tion of root and row crops, especially maize in Central 
Europe (Klik & Eitzinger 2010). 

A significant earlier beginning of silking 
(P = 0-0027) and physiological maturity (P=0-0003) 
started in 1998 and 1996, respectively, while a signi- 
ficant decrease in grain yield (P = 0-0418) commenced 
in 2002 (Fig. 5). 

Sensitivity analysis of maize production in future 
climate 

Changing the initial weather and C0 2 conditions 
of the DSSAT model, according to the results of 
different climate scenarios at the end of 21st century 
(Vucetic 2008), the sensitivity of maize growth and 
yield in the future climate was analysed. In the future 
climate, in central Croatia, global solar radiation is 
projected to increase by c. 1-07 times the present 
values (period 1949-2004), minimal temperature 
by 2 °C and maximal temperature by 4 °C, but the 
amount of precipitation is projected to decrease by 
0-92 of present values. 

Separate increases in minimal and maximal temp- 
erature and a decrease in precipitation amount, com- 
pared with present (1949-2004) climate, led to a 
prediction of maize yields being reduced by 6, 12 and 
3%, respectively (Fig. 6). A doubling of C0 2 concen- 
tration stimulated leaf assimilation, but maize yield 
was only 1% higher. Global solar radiation growth by 
7% increases evapotranspiration by 3%. By applying 
all the modifications simultaneously, with the same 
practical management and hybrids as used currently, 
the length of the maize-growing period would be 
shorter by about a month and the maize grain yield 
would be reduced by 9% compared to the long-term 
average. 

DISCUSSION 

Weather conditions and field experiment 

Despite the decrease in total arable maize area in 
Croatia, the significant increase in maize grain yield 
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Table 5. Predicted mean (MEAN), maximum (MAX) and minimum (MIN) values of particular maize 
parameters according to the DSSAT 4.0 model based on the Zagreb-Maksimir meteorological data for the period 
1949-2004. s.d. is the standard deviation and AMP L is the amplitude 



1949-2004 


Silking 
date 


Maturity 
date 


Grain yield 
(t/ha) 


Kernel 
mass (g) 


Max. LAI 

(m 2 /m 2 ) 


Biomass 
(t/ha) 


Grain N 
(kg/ha) 


Total N uptake 
(kg/ha) 


MEAN 


19 Jul 


26 Sep 


12-4 


0-46 


4-4 


22-2 


183 


257 


S.D. 


6 


18 


1-48 


0-06 


0-3 


1-5 


16 


14 


MAX 


2 Aug 


13 Nov 


14-4 


0-60 


5-0 


24-7 


218 


284 


MIN 


1 Jul 


17 Aug 


8-3 


0-30 


3-7 


17-6 


146 


232 


AMPL 


32 


88 


6-1 


0-29 


1-3 


7-1 


72 


52 
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Zagreb area (1980-2004) 
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Fig. 4. Correlation between observed and predicted silking and maturity dates for maize in the Zagreb area during the period 

1980-2004. 



(Fig. 2) is the result of sowing late-maturing hybrids 
with longer growing periods combined with improved 
agricultural management. For example, in 2008 the 
mean grain yield was c. 8 t/ha, which is the maximal 
yield during the whole period 1885-2008. According 
to recommendations by Majdak et al. (2007) autumn 
soil fertilization should be undertaken under a deep 
furrow (0-25-0-30 m) and spring fertilization should 
be carried out before sowing and reinforced after 
emergence. Thus, a grain yield of 10-15 t/ha requires 
c. 200-250 kg N/ha. In the 1999 maize field exper- 
iment, 248 kg/ha were used. As already shown, the 
weather conditions were good, especially from sowing 
to silking and the beginning of grain filling. Warm 
weather, with some drier conditions, from silking to 
maturity accelerated grain filling and physiological 
maturity, which took place around 10 days earlier. 
The DSSAT 4.0 simulation did not suggest the 
occurrence of water stress during the whole maize- 
growing period. 

However, the harvest of maize not only depends on 
weather conditions, but also on human decisions. The 
best moment for harvest of dry maize grain is when 
relative humidity is lower than 30%. In order to meet 



this condition, the 1999 harvest was rather late, on 
16 October 1999. Although the first part of September 
1999 was dry, heavy rain (82-6 mm) fell at the 
beginning of October (Fig. 3) and delayed the onset 
of the harvest. 

Thus, it is considered that the field experiment data 
were representative of the initial physiological and 
cultivation data in the study of potential crop 
production over expected climate change conditions. 
However, one of the problems in the study of the 
climate change impact on agricultural production is 
the lack of model input data, such as crop and soil 
data, for the long-term period. 

Calibration and validation of the DSSAT 4.0 

The predicted values of the DSSAT 4.0 model, 
derived from the 1999 meteorological data and 
from the pedological, physiological and cultivation 
data measured during the 1999 field experiment, are 
very similar to the observed values: beginning of 
silking and maturity, maximum LAI and biomass 
(Table 4). However, the DSSAT 4.0 simulation 
underestimated the grain yield and harvest index, 
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Fig. 5. (a, c and e) Predicted time series, significant linear trends (P<0-05) and 5-year running averages of the beginning of 
silking and maturity (days) and grain yield (kg/ha) according to the DSSAT maize model for Zagreb in the period 1949-2004. 
x is a number of years (0, 1,2,..., 55). (b, d and/) The progressive trend test for the significant linear trend in the beginning of 
silking and maize maturity (days) and grain yield (kg/ha), for the forward series u(t) (thick line) and backward series u'(t) (thin 
line) for Zagreb during the period 1949-2004. The positive u(t) points at an increasing trend, while the negative u(t) points at a 
decreasing trend. In order to identify the beginning of the possible trend, u(t) has been calculated from the first to the last 
datum, forming a progressive onward test series. The backward test series u'(t) has been formed in the same manner, 
calculating it from the last to the first term. If there is no trend, the u(t) and u'(t) curves overlap several times, whereas in the 
case of a trend, the intersection point designates the beginning of the trend, becoming significant at P<005 in case the 
absolute u(t) exceeds the ± 1-96 values (thin horizontal lines). 



and overestimated kernel mass, grain N and total N deviation of up to 0-20 of simulated yield from 
uptake. This underestimation of predicted grain yield measured yield is considered good (Alexandrov et ah 
by 0-095 is considered a good assessment because a 2001). 
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Table 6. Linear trends of predicted time series of the 
beginning of silking and maturity (days), grain yield 
(kglha), kernel mass (g), number of kernel per ear, 
maximum leaf area index (LAI max, m 2 lm 2 ), 
aboveground biomass (kglha), total nitrogen uptake 
(kglha,) nitrogen in harvest seeds (kglha), harvest 
index, evapotranspiration (mm), evaporation (mm), 
soil evaporation ( mm ) and runoff ( mm ) during the 
vegetation period according to the DSSAT maize model 
for Zagreb in the period 1949-2004. Linear trends 
significant at are in bold 



Predicted values 


Trend (/decade) 


Silking day 


-1-4 


Maturity day 


-4-5 


Grain yield (t/ha) 


-0-216 


Kernel mass (g) 


-001 


Number of kernel per ear 


-0-2 


Maximum LAI (m 2 /m 2 ) 


001 


Biomass (t/ha) 


-0122 


N grain (kg/ha) 


00 


Total N uptake (kg/ha) 


0-6 


Harvest index 


-0-01 


Evapotranspiration (mm) 


3-3 


Evaporation (mm) 


0-5 


Soil evaporation (mm) 


2-8 


Runoff (mm) 


0-6 



Predicted mean values of maize components 

The DSSAT 4.0 model offered a good assessment of 
maize phenology in the Zagreb area (Fig. 4). The 
underestimation of predicted grain yield compared 
with the actual observed yield was linked to a short, 
but very warm period during grain filling, i.e. to 
greater amplitude between the maximal and minimal 
temperature. The DSSAT model seems to indicate a 
better grain filling rate when the temperature differ- 
ence between maximal and minimal is minimized 
(Wilkens & Singh 2001). As the sensitivity range of 
these differences strongly influences the predicted 
values of grain yield, it needs to be evaluated against 
a model. 

Predicted values of phenological stages similar to 
those in Croatia, have been determined by Hunkar 
(1994) for Hungary. The beginning of the silking date 
in Hungary, as calculated by the DSSAT 2.1 model, 
takes place at the same time as in Croatia. On average, 
maturity begins about 2 months after silking, i.e. 10 
days earlier than in Croatia. The predicted values of 
grain yield and biomass for Hungary are, on average, 
smaller by 1700 kg/ha than for Croatia. However, the 
absolute value of grain yield is not so important for 
research in the climate change impact on crop 
production -the relative relationship between grain 
yield in the present and future climate is more 
important. 



Maize production in present and future climate 

Modelling of the impact of climatic changes on maize 
growth and development indicates a double C0 2 
effect. Carbon dioxide in the atmosphere directly 
affects plants by increasing photosynthesis and de- 
creasing transpiration. An increase in C0 2 in the 
atmosphere also changes the weather conditions 
affecting plant growth. Thus, an increase in air 
temperature shortens the vegetation period and 
reduces maize yield, while an increase in solar 
radiation stimulates not only leaf assimilation but 
also evapotranspiration. Taking the direct impact of 
CO2 into consideration, together with its double 
concentration in the stochastic weather generator 
Met&Roll (Dubrovsky et al. 2000), resulted in a 
9-10% increase in potential maize yield predicted for 
the Czech Republic (Zalud & Dubrovsky 2002). The 
direct and indirect C0 2 effects, combined, suggest a 
5-14% increase in potential yield. It is also shown that 
the increase in the simulated yields of other agricul- 
tural crops (e.g. winter wheat, soybean, etc.) for the 
21st century was primarily because of the beneficial 
influence of the direct C0 2 effect (Alexandrov et al. 
2002). 

It has already been noticed in Croatia that silking 
and maturity are occurring significantly earlier, and 
grain yield is significantly decreased, in the present 
climate. In previous studies (Vucetic 2006, 2008, 2009), 
it has been shown that different transient climate 
change scenarios during the 21st century, including 
only an indirect C0 2 effect, gave projections of a 
shorter growing season (34 days shorter for CSIRO- 
Mk2, 43 days for ECHAM4/OPYC3 and 44 days for 
HadCM3) and a reduction in maize yields by 10, 8 
and 15%, respectively. In this study, using both direct 
and indirect C0 2 effects in a changed climate, the 
maize-growing period has been predicted to be a 
month shorter and the maize grain yield reduced by 
9% compared with the present climate. These assess- 
ments agree with the results from neighbouring 
countries. Bacsi & Hunkar (1994) have come to the 
conclusion that the growth period of maize in western 
Hungary in changed climate would be 20^41 days 
shorter, which would result in a 7-14% smaller grain 
yield. Research in Slovenia shows that on the 
assumption of a 2 °C increase in temperature the 
cultivation area could be raised to a higher altitude 
(Kajfez-Bogataj 1993,1996,1998). It has been pre- 
dicted that the maize yield at 500 m under present 
conditions would correspond to the yield at 900 m at 
the end of 21st century. A sensitivity analysis in the 
main arable areas of the European Community (Wolf 
& van Diepen 1995) reveals that the potential yield of 
grain maize increases in the northern part of the EC, 
remains similar in the central part and decreases in the 
southern part. Drought is the main factor limiting 
the productivity of crops in southern Europe. For 
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Fig. 6. Comparison between predicted maize parameters, evapotranspiration and soil evaporation in present climate (1949- 
2004, thin horizontal lines) and in future climate (bars and thick line) for Zagreb, which is representative of central Croatia. 



example, the mean reduction of barley yield by 
terminal drought was 27% (Gonzalez et al. 2010). 
The Northern Atlantic Oscillation (NAO) exerts a 
dominant effect on air temperature in Europe and 
for example in England there has been a signi- 
ficant correlation between the winter NAO and 
summer climate since the mid-1950s (Atkinson et al. 
2008). 

This paper confirms that, in the future, Croatia 
could belong to the area of decreased maize yields. 
Some adaptation options, such as shifting to an earlier 
sowing date and selecting hybrids resistant to drought, 
could be an appropriate response to offset the negative 
effect of an increase in temperature. For example, 
using hybrids with a medium growing season would 
be beneficial for maize productivity in Bulgaria in the 
changed climate (Alexandrov & Hoogenboom 2000). 

Thus, the main conclusion from this study is that a 
direct impact of C0 2 does not indicate a considerable 
increase in grain yield in the changed climate in 
Croatia, whereas extremely warm conditions would 
have a definite influence on maize yield reduction. 
Some other crops also have a significant negative 
response to higher daily maximal temperature, for 
example, soybean (Zheng et al. 2009). 



The experience of farmers and agronomists in 
central Croatia shows that they have already been 
adapting maize production to the warmer weather 
conditions in the last decade. The usual hybrids in 
central Croatia were hybrids with a medium growing 
season but nowadays hybrids with a longer growing 
season are increasingly used. This study could help in 
optimizing and improving agricultural management 
in order to adapt for changes in climate and weather 
conditions in Croatia. 
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project COST Action 734 Impacts of climate change 
and variability on European agriculture (CLIVAGRI) 
and Climate variations and changes and response in 
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Gordon Y. Tsuji, James R. Kiniry, Josef Eitzinger, 
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